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We study the pole structure of the (-function associated to 
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in the half-line R + , subject to the singular potential gx~ 2 + 
x 2 . We show that H admits nontrivial self-adjoint extensions 
(SAE) in a given range of values of the parameter g. The 
("-functions of these operators present poles which depend on 
g and, in general, do not coincide with half an integer (they 
can even be irrational). The corresponding residues depend 
on the SAE considered. 
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I. INTRODUCTION 

In Quantum Field Theory under external conditions, 
quantities like effective actions and vacuum energies, 
which describe the influence of external fields or bound¬ 
aries on the physical system, are generically divergent 
and require a renormalization to get a physical meaning. 

In this context, a powerful and elegant regularization 
scheme is based on the use of spectral functions, such as 
the associated (-function [5j2l and heat-kernel (for recent 
reviews see, for example, f!| |(J). 

It is well known j^,^] that for an elliptic boundary 
value problem in a ^-dimensional compact manifold with 
boundary, described by a differential operator A of or¬ 
der u), with smooth coefficients and defined on a domain 
of functions subject to local boundary conditions, the (- 
function 

Us) = Tr{A~ s } ( 1 ) 

has a meromorphic extension to the complex s-plane pre¬ 
senting isolated simple poles at s = (v — j)/u ;, with 
3 = 0 , 1 , 2 ,... 

In the case of positive definite operators, the (-function 
is related, via Mellin transform, to the trace of the heat- 
kernel of the problem, Tr{e~ tA }. The pole structure of 
(a(s) determines the small-f asymptotic expansion of this 
trace [||, 

OO 

Tr{e~ tA } ~ a i ^ t U ~ u)/u , (2) 

3=0 

where the coefficients are related to the residues by 

aj(A) = Res| a _^_ J - )/w r(s) (a(s). (3) 


However, for the case of a differential operator with 
coefficients presenting singularities, less is known about 
the structure of the (-function or the heat-kernel trace 
asymptotic expansion. 

Callias |)-11| has argued that, when the coefficient in 
the zero-th order term in an elliptic, (essentially) self- 
adjoint, second order differential operator presents a sin¬ 
gularity like 1/x 2 , the heat-kernel trace asymptotic ex¬ 
pansion in terms of powers iCi-")/ 2 ( as j n (jg)) j s ill- 
defined, and an expansion including logf and perhaps 
more general powers of t ( t a with a ^ n/2) would 
be in order. In particular, considering Hamiltonians H 
with these characteristics, it has been given in a 

small-f asymptotic expansion for the diagonal element 
e~ tH (x,x) which also presents ff? - ")/ 2 logf terms, and 
where some of the coefficients are distributions with sup¬ 
port concentrated at the singularities. 


It is the aim of the present article to analyze the 
pole structure of the (-function of a Hamiltonian H de¬ 
scribing a quantum Schrcedinger particle living in the 
half-line R + , subject to a sin gula r potential given by 
V(x ) = gx~ 2 + x 2 , for a real g 12|. 

For a certain range of values of g , this Hamiltonian (a 
second order differential operator) admits nontrivial self- 
adjoint extensions! 1 in L 2 (R + ), each one describing a dif¬ 
ferent physical system. We will show that the associated 
(-function presents isolated simple poles which depend 
on g , which (in general) do not lie at s = (1 — j )/2 for 
j = 0,1,..., and can even be irrational numbers. More¬ 
over, we will find that the residues at these simple poles 
depend on the self-adjoint extension of H considered. 

This pole structure for the (-function implies a small- 
t asymptotic expansion for the heat kernel trace of the 
problem in terms of powers which (in general) are not half 
an integer. Moreover, the coefficients in this expansion 
depend on the selected self-adjoint extension. 


The structure of the paper is the following: In Section 
|n| we specify the adjoint of the Hamiltonian operator and 
in Section fin we determine its deficiency subspaces. The 
Hamiltonian self-adjoint extensions are characterized in 
Section |v|, and in Section [v| is described the correspond- 


1 The existence of nontrivial SAE for this kind of singular 
potential has been pointed out in j^| . SAE with more general 
singular potentials have also been considered in ]2^,p3|. 
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ing spectrum. In Section VI we give an integral represen¬ 
tation for the ^-function of each SAE of the Hamiltonian 
and in Section VII w e disc uss the structure of its sin¬ 
gularities. In Section VIII we analyze some particular 
cases, and we establish our conclusions in Section IX. 
Appendix |a| is devoted to the construction of the clo¬ 
sure of H, and in Appendix |b] we outline the necessary 
asymptotic expansions. 


II. THE HAMILTONIAN AND ITS ADJOINT 

Let us consider the operator 


H = ~h +V ^ 

( 4 ) 

V(x) = zL+x 2 , 

X A 

( 5 ) 


densely defined on the domain T>(H) = Cq°(R + ), the lin¬ 
ear space of functions ip(x) with continuous derivatives of 
all order and compact support non containing the origin. 
It is easily seen that H is a symmetric operator. 

In order to construct the SAE |2(J of H we must get 
its adjoint, H ', and determine the deficiency subspaces. 

The operator H 1 is defined on the subspace of square- 
integrable functions ip(x) for which (’ll), Hip) is a con¬ 
tinuous linear functional of <p £ T>(H). This requires 
the existence of x( x ) £ L 2 (R + ) such that (ip, Hip) = 
(x,<p)y ip £ 'D(H). If this is the case, then x( x ) is 
uniquely defined, since V(H) is dense in L 2 (R + ) and, 
by definition, H^ip = %. 

For ip £ and \/ ip £ H(H) we have 

(ip, Hip) = f 0 °° ip(x)*(-<p"(x) + V(x)ip(x))dx = 

( 6 ) 

= {{-V + v(x)ip),ip) = ( x, p)> 

where the derivatives of ip are taken in the sense of dis¬ 
tributions. 

Equation (|b|) implies that ip"(x) = V(x)ip(x) — x( x )> 
is a locally integrable function. Then, its primitive ip'(x) 
is absolutely continuous for x > 0. 

Therefore [Q, the domain of H 1 is the subspace of 
square integrable functions having an absolutely contin¬ 
uous first derivative and such that 


III. DEFICIENCY SUBSPACES OF H 

To compute the deficiency indices ||(J of H, n± = 
dim/C±, we must solve the eigenvalue problem 

H^(p x = ~<Px(x)+V(x)<px(x) = X(p\, (8) 

for (px £ V(H') and A £ C, with the imaginary part 
3(A) ^0. 

By means of the following Ansatz (suggested by the 
expected behavior of the solutions of (g) for x —> 0 + and 
x —> oo), 

(p = x a e~^F(x 2 ), (9) 

with 

a = 1/2 + x/ g + 1/4, (10) 

we get from (|j) the Rummer’s equation for F(z): 

zF"(z) + (b - z)F'(z) - aF(z) = 0, (11) 

where a = (2a + 1 — A)/4 and b = a + 1/2. 

For real a, we have g > —1/4 and a > 1/2. In this 
case it can be seen Q that the only solution of eq. 
leading to a square-integrable at infinity solution of eq. 
(||) is given by the Kummer function F(z) = U(a;b;z). 
Then, the eigenfunctions of H 1 are proportional to 

(px(x) = x a e~^~ U ( ^ a —-^a + ^x 2 ^. (12) 

We must now study the behavior of <px near the origin, 
where U(a;b;z) behaves as z~ a ( 1 + 0(1/z)) |25|| . We 
must consider two different regions for the parameter a. 

For a > 3/2, <p x £ L 2 (R+) o a = (2a + 1 - A)/4 = 
—n, with n £ N. As a consequence, if A ^ R, (px ^ 
L 2 (R + ), and the deficiency subspaces are trivial. 

This means that, for a > 3/2, H is essentially self- 
adjoint, and its discrete spectrum is given by the condi¬ 
tion —a £ N, i. e. 

A n = 4n + 2a + 1, (13) 

with n = 0,1,2,... The corresponding eigenfunctions are 

(p n = x a e~~U (~n;a+ ^;x 2 j . (14) 


H*ip(x) = — ip"(x) + V(x)ip(x) £ L 2 (R+) (7) 

(without requiring any boundary condition at x = 0). 

In the next Section we will determine the deficiency 
subspaces of H, K,± = Ker(IU i). 


On the other hand, for 1/2 < a < 3/2, one can see Q 
that (px £ T 2 (R + ),VA £ C. Then, the deficiency sub¬ 
spaces IC± are one-dimensional, and the deficiency indices 
n± = In this region, H admits different self-adjoint 
extensions. 


2 This is in accordance to Weyl’s criterion [G0| according to 
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IV. SELF-ADJOINT EXTENSIONS OF H 


Since n+ = 1 = n_ for 1/2 < a < 3/2, there exists 
a one-parameter family of self-adjoint extensions of if, 
which are in a one-to-one relationship with the isometries 
from K.+ onto /C_ f 20 j . 

The deficiency subspaces /C+ and /C_ are generated by 
(j) + = 4>\-i and 4>- = = </+, respectively. Then, 

each isometry if 7 : /C+ —> /C_ can be identified with the 
parameter 7 £ [0, 7 r) defined by 


« 7 <£+ =e " 2i7 <t>-. (15) 

The corresponding self-adjoint operator, if 7 , is defined 
on a dense subspace [|o| 

P(if 7 ) C 2?(ff t ) =P(if)©/C+ ®AC_, (16) 

where if is the closure of if. Functions £ T>(H 7 ) can 
be written as 


— </>o + A (/>+ + e , (17) 

with <j>Q £ 1 ?(if) and A a constant. Since if 7 is a restric¬ 
tion of ift, we have 


H 1 (j> = H^<j> = H ( j) 0 + iA((t> + ~e- 2i ' y (t>-) . (18) 

In the following we take g>0=>l<a<3/2. As we 
will see, condition (pd]) determines the behavior of <f> £ 
1 ?(ff 7 ) near the origin. Taking the logarithmic derivative 
of <j> we get, 


V e*</> t 0 +2AK(e* <(>'+) 

c\> e i7 </ 0 + 2A3?(e i 'i'<(> + )' ' ’ 

In this expression, the terms coming from <p + give the 
leading contributions for small x. In fact, in Appendix 
|A| we show that <j) 0 (x) = o(x a ) and 4>' n (x) = o(x“ _1 ). 
Then, for the right hand side of eq. ([H]) we get |2f| (see 
eq- ( 0 ), 


which, for continuous V(x), H is essentially self-adjoint if and 
only if it is in the limit point case, both at infinity and at the 
origin. 

In addition, if V(x) > M > 0, for x large enough, then H 
is in the limit point case at infinity. In consequence, in the 
present case H is essentially self-adjoint if and only if it is in 
the limit point case at zero. 

In particular, for positive V(x) (g > 0), if V(x) > 3/4 x~ 2 
for x sufficiently close to zero then H is in the limit point case 
at the origin. On the contrary, if V(x) < (3/4 — e) x~ 2 , for 
some e > 0, then H is in the limit circle case at zero. 

This confirms our results concerning the self-adjointness of 
H in the different regions of the parameter g. 


&’(x) 

</>(x) 


1 — a 
x 


+ ( 2 a 


V(a —i) 


cos ( 7 - 71 ) 
cos ( 7 - 72 ) 


2a—2 


+ o(x Ja ~ 2 ) 


( 20 ) 


where we have called 71 = arg{T[(— 2a + 3 — *)/4]} and 
72 = arg {T[( 2 a + 1 - i)/4j}. 

Thus, the limit of eq. (jl^) for x —> 0 + gives the appro¬ 
priate boundary condition for the functions in the domain 
of the particular SAE . As we will see, this boundary con¬ 
dition will finally determine a discrete spectrum for H^. 


V. THE SPECTRUM 


The boundary condition specified in eq. (|2^) character¬ 
izes the domain of a particular SAE of the operator H , 
H 1 . In order to determine its spectrum, we must find the 
solutions of (||), (j)\ as given in ( 0 ) with A £ R, which 
satisfy this boundary condition. Their behavior near the 
origin is given by (see eq. (0 and p5||), 


^aO) 

(t>\ (x) 


1 — a T(s — a) 

-+ (2a - l )-/2 - 7 — x 


p |~ —2a+3—A j 


T(a-i)' 


: + o(x 2a ~ 2 ). 


( 21 ) 


Comparison of eqs. 
to 


and ( 21 ) immediately leads us 


—^- *L-=/3( 7 , k ), ( 22 ) 

where we have defined the parameters 

k = —-— = j (2 + y/l + 4g) £ [3/4,1) 

4 4V J (23) 

/3(7, k) = cos (7 - 71)/ cos (7 - 72). 

Eq. (^ 2]) determines a discrete spectrum for each SAE. 
In Figure 1. we plot both sides of eq. (|22|) as a function 
of A, for k = 4/5 and (3=1. The abscissae of the in¬ 
tersections of this two functions give the corresponding 
spectrum. 

Notice that each SAE can equivalently be character¬ 
ized by (3 £ R U {—00}. Then, we will also use the 
notation to design this SAE. 
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(28) 



FIG. 1. F( A) = —=—as a function of A, for re = 4/5. 

r (, 1 — K ~ T ) 

The solutions of F( A) = 0 give the spectrum of the SAE 
identified by 0. 

The spectrum of H(p) is bounded from below, and 
presents a negative eigenvalue for those SAE character¬ 
ized by (3 > r(«)/r(l — re) (even though the potential 
V(x) > 2 yfg > 0). Moreover, there is no common lower 
bound; instead, any negative real is in the spectrum of 
some SAE. 

For any value of g , there are two particular SAE for 
which the spectrum can be easily worked out (see eq. 

• For j3 = 0 the spectrum is given by 

A„ = 4(n + 1 - re), (24) 

with n = 0,1, 2 ,... 

• For (3 = — oo the spectrum is given by 

A„ = 4(n + re), (25) 

with n = 0,1, 2 ,... 

For other values of (3, the eigenvalues grow linearly 
with n, 


4 (n — 1 + re) < X n < 4 (n + re). (26) 


The case with g = 0 

It is instructive to consider the particularly simple case of 
the harmonic oscillator in the half-line, for which there is no 
singularity in the potential. Indeed, for g = 0 (a = 1 or 
re = 3/4), the boundary condition (eq. (|2c|)) reads, 

■TT-y = —2/3 + 0(x) (27) 

<P{x) 

or, equivalently, 


liin + 2f3cf>(x)\ = 0, 

which corresponds to Robin boundary conditions at the ori¬ 
gin. Dirichlet and Neumann boundary conditions are ob¬ 
tained for /3 = —oo and 0 = 0, respectively. 

Let’s now study the eigenfunctions and eigenvalues of the 
self-adjoint extensions of H corresponding to different values 
of 0. 


Dirichlet boundary conditions (0 = —oo) 

Since re = 3/4, the eigenvalues (see eq. (ji^)) are given by 
A n — 4 n + 3, (29) 

where n = 0,1, 2,... 

Since the Hamiltonian (eq. (f)) corresponds in this case 
to a particle with mass m = 1/2 and frequency u> = 2, the 
eigenvalues of this SAE can be written as X n = ui[{2n + 1) + 
1/2], which coincides with the spectrum of the odd parity 
eigenvectors of the harmonic oscillator on the complete real 
line. 

In fact, the eigenfunctions are given by (see eq. (^2j) and 

II), 

= 2- 2n - 1 e-^-H 2n +i{x). (30) 


Neumann boundary conditions (0 = 0) 

In this case, the eigenvalues (see eq. (^)) are given by 

A„ = 4n + 1, (31) 

where n = 0,1,2,.... This eigenvalues can be written as 
A n = w(2n + l/2), which coincides with the even parity sector 
of the harmonic oscillator spectrum on the complete real line. 
The eigenfunctions are now given by (see eq. (|l^)), 

0 n = 2- 2n e-^H 2n (.x). (32) 


Robin boundary conditions (0 ^ 0, — oo) 

For finite 0^0, the eigenfunctions are given by (eq. (|l^)), 

<t, x =xe-^u(^^-,^,x 2 y (33) 

and the corresponding eigenvalues are determined by the 
trascendental equation 

F(3zA) 

pfey =0(7,3/4). (34) 

Notice that, for general Robin boundary condition, the 
ground state is negative (less than the minimum of the po¬ 
tential) for 0 > r[|]/F[i], 
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VI. THE INTEGRAL REPRESENTATION FOR 
THE C-FUNCTION 


The spectrum of each SAE of the operator H in ( 0 ) 
is determined by eq. (f22|), for any given (3 £ [— 00 , 00 ). 
In this section, we will study the pole structure of the 
associated ^-function, defined as 

= (35) 


Notice that, since the eigenvalues grow linearly with n 
(see eq. (pbj)), C/?( s ) is analytic in the half-plane §i(s) > 1. 
For finite /3, let us define the holomorphic function, 


/(A) 


1 

r(i-K-S) 


(3 

r («-*)’ 


(36) 


with | < k < 1. The eigenvalues of the self-adjoint 
operator H^ correspond to the zeroes of /(A) which, 
consequently, are all real. They are also positive, with 
the only possible exception of the first one, according to 
the discussion in the previous Section. 

Moreover, the zeroes of /(A) are simple. To prove this, 
let’s assume the converse is true, i. e. there is a A £ R 
such that /(A) = /'(A) = 0. Taking into account that 


1 

4 


/'(A) = 

[V»(l-K- 


_ ^(l-«~|) 

' 4 r(i- K -|) 

IHfci)1 


r(i-K-l) 


-0 


V , ( k -t) _ 
4 r( K -A) 


+ ip(n- A/4)/(A)| , 


(37) 


where ip{ z ) is the polygamma function, we see that our 
assumption requires that 


ip (1 — k — A/4) = ip (k — A/4), (38) 

which is not the case for any A £ R, if f < « < 1. 

Therefore, the ([-function can be represented as the 
integral on the complex plane 

where C is a curve which encircles the positive zeroes of 
/(A) counterclockwise. In eq. (|39|) , Q(y) = 1 for y > 0 
and 0(y) = 0 for y < 0. 

Let us consider the dominant asymptotic behavior of 
the quotient 


/'(A) 

/(A) 

[*(i-.-i)-*(.-i a + i («> 


4 1 -D 


r(i-tt-f) 

r («-T) 


For | arg(—A)| < 7r and |A| —> 00 , it is sufficient to write 


V>(k- A/4) = log(-A) + 0(1), (41) 


ip{l — k — A/4) — ip(K — A/4) = 0(A _1 ), (42) 


r («-|) 


0(A 1_2k ). 


(43) 


Consequently, for 5ft(s) > 1 the path of integration in 
(p9|) can be deformed to a vertical line, to get 




zoo-j-0 f'(\\ 

^T^dX + his), < 44 ) 

-ioo+0 J{ A ) 


where h(s) (the contribution of the negative eigenvalue, 
if any) is a holomorphic function. 


VII. POLE STRUCTURE OF THE C-FUNCTION 


The integral in eq. ( fl4|) defines C/ 3 (s) as an analytic 
function in the half-plane 9?(s) > 1, which can be mero- 
morphically extended to the whole complex s-plane. It 
can be written as 

. 1 f i0 °f'W ,-S „ 

C/3(s) - 2** l /(A) 


1 

n 


A _s dX + hi(s) = 


J-i 00 /(A) 

e -isn/2 roo ft ( ifi) 

2 ?r A f{in) 


(45) 


y s dy— 


e”/ 2 f'{—iy) _ SJ , , . 

' y dy + hi(s), 


2tt 


f(~¥) 


where hi(s) is a holomorphic function. 

In Appendix [b] we work out the asymptotic expansion 
of /'(A)//(A), which is given by 

f'(\) 1 1 00 

J[ ’ -log(-A) + -^c fc ( K )(-A)- fc + 


/(A) 


k=0 


(46) 


(- A ) 


— N(2K — l) — 2n—l 


N—l n=0 


where the coefficients Cfc(ft) are polynomials in n whose 
explicit form is not needed for our purposes, and 


= — (4 2k-1 /3) JV (^2k — 1 + b n (K,N), 
with b n (n,N) defined in eq. (|B7|) (see also eq. (|B8|)). 


(47) 
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As can be seen from eq. (|46|), the asymptotic expansion 
of /'(A)//(A) contains the logarithmic term |log(—A), 
and a series of non positive integer powers of A, both 
coming from the /^-function in the last term in the right 
hand side of (|4(]). There is also a series of decreasing re- 
dependent powers of A, which comes from the first term 
in the right hand side of (fft|. 

For the dominant logarithmic term we get from ( phi) 
log (e~ 1 ^ p) + e 1 ^ log (e l %p) ] p~ s dp 



sin (Ar) cos(^) 

8 (s — 1) 4-7T (s — l) 2 


1 1 

4 (s — 1 ) 


h 2 {s), 


(48) 


with a /3-dependent residue given by 


R- eS (')ls=-XV(2K-l)-2n — 
(~1) N 

- Cn,ti(kiP) sin(27rAfre). 


( 52 ) 


This is our main result, establishing the existence of re¬ 
dependent poles of the ((-function which, in general, are 
not located at a half an integer value of s. Moreover, the 
residues depend on the SAE considered. 

Finally, notice that when re is a rational number, there 
can be several (but a finite number of) pairs ( N , n ) con¬ 
tributing to the same pole. They must satisfy 


n — n! 
N-N' 


«=— (-1/2,-1/4], 

z q 


(53) 


where h 2 (s) is holomorphic. The analytic extension of 
this term presents a unique simple pole at s = 1, with a 
residue equal to 1/4. 

The remaining terms in the asymptotic expansion of 
/'(A)//(A) are of the form A,-(—A)~ J , for some j > 0 
(see eq. (|b[)). Replacing this in eq. (|4E|) we get 


A 

27T 


,-if (s-j) gif (s-j) 


dp = 



(49) 


_ Aj sin(7r j) 1 

— T-, rr + n-3{s), 

n s-(l-j) 

where h^{s) is holomorphic. 

So, from each power dependent term in the asymptotic 
expansion of /'(A)//(A), proportional to (—A) _J , we get 
a unique simple pole at s = 1 — j, with a residue given 
by “(A? AO sin(7rj). 

Notice that this residue vanishes for integer values of 
j. In particular, this is the case for all the contribution 
coming from the asymptotic expansion of _0(re — A/4) in 
the last term in the right hand side of eq. ( po| ) , except for 
the first one, the logarithmic term leading to eq. (p8|). In 
fact, this is the only singularity present in the (3 = —oo 
and (3 = 0 cases (see eq. ©)■ 

But in general, for | < re < 1, there are also poles at 
non integer values of s, as follows from (|4fj). 


In conclusion, besides the pole at s = 1 with residue 
1 /4, for each pair of integers 


where p, q £ N. 

On the contrary, when re is irrational the poles com¬ 
ing from different pairs ( N,n ), also irrational, are not 
coincident. 


A. Poles and residues of </g(s) 

Let us recall that the logarithmic term in the expansion 
( fl6| ) leads to a pole at s = 1 (see eq. (|4^)) with a residue 
given by 


Res(-)| S=1 = i, (54) 

independently of the SAE considered. 

The other poles can be organized in sequences char¬ 
acterized by the integer N = 1,2,... In each sequence, 
successive poles differ by —2. 

For example, the poles corresponding to the pairs (N = 
1, n), with n = 0,1, 2,..., are located at (see eq. (|5l|)) 

1 , , 

— 1 — 2n<s = l — 2re — 2n< — — — 2 n, (55) 

and have residues given by 

Res (■)I s=i— 2 K— 2 n = - Cl ’ n ( K ’® sin(27rre). (56) 

7T 

Similarly, the poles arising from the (N = 2, n > 0) 
terms in the asymptotic expansion (|f|) are located at 

— 2 — 2n < s = 2 — 4re — 2n< —1 — 2 n, (57) 


(N, n), with N = 1,2,3,..., and n = 0,1,2,..., (50) and have residues given by 


the /'-function of the SAE of H characterized by the pa¬ 
rameter /3, //j(s), has a contribution with a simple pole 
at the negative value 

N 

s = —_/V(2re - 1) - 2n e (-N - 2n ,- 2 n\, (51) 


R eS (")ls=2—4re—2n = S in(47TK). (58) 

7T 

Notice that the poles in the N -th sequence have 
residues proportional to [3 N (see eq. (|T7|)) . 
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Finally, let us stress that a pole of ^g(s) at a non integer 
s = — N(2n — 1) — 2 n, as in (|5l|), implies the presence of a 
term in the small-f asymptotic expansion of Tr{e -tff<3) } 
of the form 


a ±N (2n-l)-\-2n 

A [N(2K,-l)+2n\ t > 

with a coefficient related to the residue by 


A 


[iV(2K-l)+2n] 


= r(-JV(2/« - 1) - 2n) Res C/3 (s) | s =_jv (2 «- 1) - 2 r l ■ 


(59) 


(60) 


B. (-function singularities from the asymptotic 
expansion of the eigenvalues 

The singular behavior found for ^g(s) can be confirmed 
(at least for the first few poles) by determining from ( f^ ) 
the asymptotic expansion of the eigenvalues for n 3> 
1. Indeed, one can make the Ansatz 


VIII. PARTICULAR CASES 

In this Section we will show how our results reduce to 
the usual ones for g = 0 (when there is no singularity in 
the potential). We will also show that, for /? = 0 and 
(3 = — oo, the ((-function presents a unique simple pole. 


A. The (3 = 0 and (3 = —oo SAE 

The ((-function for the SAE characterized by j3 = 0 and 
/? = — oo can be exactly evaluated, since in these cases 
the spectrum was explicitly computed in eqs. (^4j) and 
( p5| ) respectively. We get 

OO 

Co00 = 4 _s ©(n + 1 - k)~ s = 4“ 5 C(s, 1 - k), 

n—0 

(64) 

OO 

C-oo(s) = 4 -s ^2(n + k)~ s = 4~ s C(s, k), 

n =0 


= \ -K + n + e, (61) 

and self-consistently determine e through successive cor¬ 
rections. For the first terms we get 

13,71 ~ 1 — k + n + — sin(2 7r k) n 1 ” 2 K + 

4 7r 

+— (l — 3 k + 2 k 2 ) sin(2 7r/;)n _2K - (62) 

7T 

sin(4 7TK) n 2_4ft + ..., 

2 7r 


where we have retained only powers of n g reater than 
—2. This leads, for the ((-function in eq. (35), to 


C/3(s)~4 s C(s)+s4 s (k- 1)C(s + 1) + 


+s (s + 1) 4 


- S ( K ~ 1) 


C(s + 2)— 


where C( s i?) is R ie Hurwitz ((-function, whose analytic 
extension presents only a simple pole at s = 1, with a 
residue Res((s, g)| s =i = 1. This leads, in both cases, to 
a unique simple pole for the ((-function at s = 1, with a 
residue equal to 1/4, in agreement with eq. (|48|). 

In fact, from eqs. (|5^) and (^) it is evident that all the 
residues corresponding to negative poles vanish for /3 = 0. 
On the other hand, for (3 = —oo, /'(A)//(A) reduces to 
A/4) (seeeq. ©) , and the only term leading to a 
singularity is the logarithm in the asymptotic expansion 
(fiff), as already discussed (see eq. ©)). 


B. The harmonic oscillator in the half-line 

For the harmonic oscillator in the half-line (g = 0 or 
k = 3/4) we still find a simple pole at s = 1, with residue 
1 /4 (the only singularity for Dirchlet or Neumann bound¬ 
ary conditions, as previously discussed). 

For finite (3, the remaining singularities are located at 
(see eq. ©), 


-s4 s — sin (2 n 7r) £(s + 2 k) — 
7r 


N 

S = -y - 2n, N = 1,2,3,..., n = 0,1,2,..., (65) 


—s (s + 2 k) 4 s — (k — 1) sin(2 7r k) C(1 + s + 2 k) + 


(3 2 

+s 4 s -— sin(4 7r«) £(s — 1 + 4 k) + ..., 

2 7r 


(63) 


with residues given by (see eq. (|52|)), 

Res(-)| s= _f_ 2 „ = 

= ^ ^ C N ,n (« = 3/4, /3) sin > 


( 66 ) 


where ^(z) is the Riemann ((-function, which presents a 
unique simple pole at z = 1, with a residue equal to 1. 

This result shows a pole structure in agreement with 
the one previously described. 


which vanish for even N. 

Then, each pole (except for the first one, at s = 1) 
corresponds to a negative half-integer, 
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s = —k — 1 / 2 , fc = 0 , 1 , 2 .... 


(67) 

Moreover, it is clear that for finitely many pairs ( N,n ) 
satisfying N + An = 2k + 1, the corresponding poles lie 
at the same point. 

Therefore, the residue of /g(s) at s = —k — 1/2 must 
be computed by adding all these contributions, charac¬ 
terized by N = 2 (k — 2 n) + 1, with n = 0,1, 2,..., [k/2]. 
We get 


Res(C/3(s))| s= _ fe _i = 

/•nk+i [ fc / 2 ] ( 68 ) 

= -- C[2(fc-2n) + l],n ( K = 3/4,/?). 

n—0 

For example, for k = 0, the residue is 
Res (C/3( s ))] s =-A =--Ci,o(k = 3/4,0) = - (69) 

27r 7 r 

and, for k = 1, 

Res(C/3(s))l s -_3 =-C' 3i0 (k = 3/4,/?) = -- /3 3 . (70) 

2 n 7r 


IX. CONCLUSIONS 

In this article we have analyzed the pole structure of 
the (/function of the Hamiltonian describing a quantum 
Schroedinger particle living in the half-line R + , subject 
to the singular potential V(x) = gx~ 2 + x 2 . 

We have specified the domain of the adjoint of the 
Hamiltonian, H '. and determined the deficiency sub¬ 
spaces of H , initially defined on Cg°(R + ). We have 
shown that, for —1/4 < g < 3/4, H admits nontriv¬ 
ial self-adjoint extensions which depend on a continuous 
real parameter 0. 

For computational convenience, we have limited our 
analysis to the range 0 < g < 3/4. 

Once determined the closure of H (studied in Ap¬ 
pendix 0 ). we were able to characterize each SAE 3) 
by the behavior (singular, in general) near the origin of 
the functions in the corresponding domain of definition. 
This relation also allowed for the identification of the 
spectrum of 3) with the zeroes of an analytic function 

/(A). 

The asymptotic expansion of /(A) (outlined in Ap¬ 
pendix |b|) led to the determination of the poles and 
residues of the /'-function associated with H^y //?(s). 

We have shown that the poles of /g(s) can be organized 
in sequences characterized by an integer N = 1,2, 3,..., 
and are located at s = —N(2 k — 1) — 2n, with n = 
0,1, 2,... and k = (1 + g + l/4)/2 e [3/4,1). Notice 
that these values of s are not, in general, half an integer 
(which are the expected positions of the poles for a sec¬ 
ond order differential operator with smooth coefficients 


on a compact segment), and they are irrational numbers 
for irrational values of k. 

We have also found that the residues depend on the 
parameter 0 characterizing the SAE H^y 

We have confirmed this /-function pole structure (for 
the first poles) through the comparison with the results 
obtained from the asymptotic behavior of the eigenval¬ 
ues. 

These results also imply that the small-t asymptotic 
expansion of the heat kernel of 3) contains powers of 
t which (in general) are not half an integer, and that the 
corresponding coefficients depend on the SAE. 

Finally, several particular cases were analyzed, finding 
that our results are consistent with the known ones. In 
particular, for the harmonic oscillator (g = 0) in the half 
line, subject to any local boundary condition at x = 0, 
the poles lie at half-integer values of s. 


A final remark is in order: Notice that the unusual pole 
structure previously described is a consequence of having 
a potential with a moderate singular behavior near the 
origin. In fact, for g > 3/4, where the Hamiltonian H is 
essentially self-adjoint due to a stronger singular behavior 
of V (ad, the /-function simply reduces to 4 _s /(s, k) (see 
eq. (|l3|)), which presents a unique pole at s = 1 with 
residue 1/4. 


A similar pole structure is obtained for the Hamilto¬ 
nian /-function of charged Dirac particles living in (2+1)- 
dimensions, in the presence of both a uniform magnetic 
field and a singular magnetic tube with a non-integer 
flux. This problem was considered in |26|, where it was 
shown that the Hamiltonian restricted to a critical angu¬ 
lar momentum subspace admits nontrivial SAE, whose 
spectra are determined by a trascendental equation sim¬ 
ilar to (p2|). These results will be reported elsewhere. 
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APPENDIX A: CLOSURE OF H 

In this Section we will justify to disregard the contri¬ 
butions from the functions in the domain of the closure 
H to the boundary condition, eq. (pii|). Indeed, we will 
show that if (j> £ 77(77) then 

(j)(x) = o(x a ) and <fi'(x) = o(x“ -1 ) (Al) 

near the origin, for any a < 3/2. 





In order to determine the closure of the Hamilto¬ 
nian’s graph we must consider those Cauchy sequences 
in V(H) = C(f (R+), {</? n }„ G N, such that {Hp n } ne N 
are also Cauchy sequences. Notice that, since the coef¬ 
ficients in H are real (see eq. ®)) , we can consider real 
functions. 

Let us call p = p n — p m , with n, m £ N. Then ip —> 0 
and Hp —> 0 as n, m —> oo. 

Consider first the scalar product 

pOO . 

Hip) = J if y—ip" + ^2 dx = 

p oo 

= / (V 2 + -4 ^ + x V ) dx - ii^ii \\ H v\\ ->■ o 

Jo v x ' 

(A2) 

for n, m —> oo. Therefore, for g > 0, we conclude that 
Pn(x) 


Wni X )} 


nEN, 


and {a;(/?„(x)} ne N (A3) 


nEN 


are also Cauchy sequences. 

We will now prove the following 
Lemma: Let {cp n }neN be a Cauchy sequence in T>(H) = 
C“(R + ) such that, for g > 0, 1 < a < 2 and g ^ (a 2 — 

l)/4, 

\J, 


Pn{x) 


, and 


nEN 


<Pn{x) 

ry>a 1 


ne N 


are also Cauchy sequences. Then, 


Pnjx) 

rf.l+a/2 


and 


nEN 


<Pni x ) 

X a / 2 


(A4) 


(AS) 


nEN 


are Cauchy sequences too. 

Proof: As before, let p = p n — p m . First notice that, 
for 1 < a < 2, 


(A6) 


^a; 1 a / 2 p(x)) dx < J [p{x)) 2 dx+ 

/ OO 

(xp(x)) 2 dx < ||<p(:r)|| 2 + ||£cV?(a:)|| 2 - 

Then, from (f\3l), we see that \ x l ~ a l 2 p n {x) ? is also 

- I J nEN 

a Cauchy sequence. 

A straightforward calculation shows that 


<,<o+1,u ^) 2 +(- i -“ /2 ^)) 2 U- 


(A7) 


Similarly, 




a — 1 


p'{x) 


CL -hi 


V(x) 


a — 3 


X a /2 

) + 

N 

\ 2 ] 

p(x)^ 

) > dx 


(A8) 


Now, taking into account that the sum of fundamental 
sequences is also a Cauchy sequence, we see that 

+ («) 
when 77., m —> oo , for any pair of real numbers A and B. 


The coefficients of 


^MV and ( ^ X) 

x a/2 J \ X 1+a / 2 

expression of the scalar product in (|A9|) are 


in the 


A-B 


a — 1 


, . a(a + 1) . „ 

A g- o + B 9 


and 

CL A 1 


(A10) 


respectively. One can see that, by choosing one of them 
as zero, the other is non vanishing (except for g = (a 2 — 

l)/4). 

It is easily seen that these results prove the Lemma. 


For o = l, from (A3) and the Lemma, we conclude 
that 


<Pn(x) 

x 3/2 


and 


nEN 


¥>n( x ) 

x l/2 


(All) 


nEN 


are Cauchy sequences. 

Let us first suppose that g is an irrationa l num ber. 
Then, applying iteratively the Lemma from (All) one 
can show that, for any positive integer k, 


<Pn(x) 


2[l-(l/2) fc ] 


and 


nEN 


<Pn( x ) 


2[l-(l/2) fc ]-l 


(A12) 


nEN 


are Cauchy sequences. 

Finally, for any given e > 0 there are integers k\ and 
&2 such that (l/2) fcl < e < (l/2) fc2 . Taking into account 
that 


1 


< 


1 


2_£ ~ x 2[l-(l/2) fe i] 


< 


2_£ ~~ 2[!-(l/2) fe2 ] 


for 0 < x < 1, 


for x > 1, 


Pn{x) 

r 2 — s 


(A13) 


is a 
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one immediately concludes that 
Cauchy sequence. 


nEN 





































A similar conclusion is easily obtained for 

Vn{ x ) 


r l— £ 


n£N 


Let us n ow su ppose that g is a rational number. Then, 
from (A3) and (All) it is seen that we can choose an 
irrational a £ (1,3/2) from which the Lemma can also 
be applied iteratively to arrive to the same conclusions. 

In the following we will consider the behavior of the 
functions near the origin. 

For any e > 0, we can write 


X®) = [ (y a <p{y))' dy = 

Jo 




(A14) 


So, for x < 1, a < 3/2 and e small enough, we have 


A 2 /, 

<p(y) 

r 1 

y 2-e 


v'{y) 


,1—e 


^n.m —»oo 


0. 


(A15) 


Therefore, the sequence {x~ a i/? n (*)}neN, with a < 3/2, 
is uniformly convergent in [0,1], and its limit is a contin¬ 
uous function vanishing at the origin, 


lim (x a ip n {x)) = x a 


lim (x “ 0(x)) = 0. 

x— > 0 + 


(A16) 


(A17) 


In particular, for a = 0 we have the uniform limit 

lim tp n (x) = <j>(x), (A18) 

n—>oc 

which coincides with the limit of this sequence in 
L 2 (R+). 

On the other hand, we can also write 

[ y~ a+1 H<p(y) dy = -x~ a+1 tp'(x)+ 

Jo 


y-“ +t -’u-a- 


(A19) 

y y ) 


+ f y a+2 yv(y)dy. 
Jo 


Therefore, for x < 1, a < 3/2 and e sufficiently small, we 
have 


|ar a+ V(a;)| < (j y 2( ~ a+1) dy^j 

<p’(y) 


1/2 


1/2 


,2(-a+l-e) 


dy J 


a — 1| 
1/2 


il — £ 


\\ H( p(y)\\- 
<p(y) 


r 2 —e 


+ y2{ ~ a+2) dy ) 


7 n,m —>-oo 


o. 


(A20) 


Consequently, the sequence {x _Q+1 with 

a < 3/2, is uniformly convergent in [0,1], and its limit is 
a continuous function vanishing at the origin, which we 
write as x~ a+1 xi x ) : 

lim (x~ a+1 ip' n (x)) =x~ a+1 x(x), (A21) 


lim (x a+1 x ( x )) = 0. (A22) 

x— >0+ 

In particular, for a = 1 we have the uniform limit 

lim <Pn{x) = X( x ), (A23) 

n —>oo 

which coincides with the limit of this sequence in L 2 (R + ) 
(see (©). 

Let us now show that x( x ) = Indeed, for x < 1, 

we have 


<t>i x ) - / x(y) dy 

Jo 


< 


< - Pn( x ) | + 


(x(y) - v’niy)) dy 


< 


(A24) 


< \H X ) - <Pn( x ) | + ||X - <Pn II *n —>oo 0. 


So, 4>{x) is a differentiable function whose first derivative 

is x( x )- 


Equations (A17) and (A22) imply that, given ei > 0 
and a < 3/2, 


|</>(x)| < £1 x“ and |</>'(x)| < eq x° 


(A25) 


if x < 6, for some 6 > 0 small enough. This proves our 
assertion. 


APPENDIX B: ASYMPTOTIC EXPANSIONS 

In this appendix we will compute the asymptotic ex¬ 
pansion for /'(A)//(A) as given in eq. ©• 

The asymptotic expansion for the polygamma function 
appearing in the right hand side of eq. ([hj) can be easily 
obtained from Stirling’s formula pB||, 
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V>(k. - A/4) - log (—A) +5 >(k)(-A) 


—k 


(Bl) 


i=0 


where the coefficients Ci(k) are polynomials in k which 
we will not need to explicitly know for our purposes. 

On the other hand, taking into account (|43|), we can 
write asymptotically for the first term in the right hand 
side of eq. ( |40| ) 

[V>( 1-K - j) £)] 


1-/3 


r (—5) 


E oo /qAT 
7V=0 ^ 


r(.-A) 


-i AT 


ty(i-«-$)-^(*-|)] = 


E oo oN 
N =0 ^ 


r («-D 


AT 


4 +6q lo s 


F («-t) 


_ 'spoo 4 d 

Z^N=0 N d{- A) 


r (—4) 


AT 


(B2) 


From the Stirling’s formula p5fl we get 


log 


r(*-£) 

{ OO 

a m(«)(-A) 

m =1 


(1 — 2k) log(--) + 


(B3) 


-2m 


where the coefficients in the series are given by 

m 

x [(1 - K ) 2m + - 2 ” 1 ] + (2m + 1) £ x 


2m — 1 

2p-2 


i 2(m-p)+l _ ^ _ K j2(m-p) + l 


(B4) 


Then, 




l N 


r («-£) 

— N(2k— 1) oo 


(B5) 


( a \ — iV — I ) OO 

- 4 ) ^6 n («,7V)(-A)- 2 ", 

' n =0 


00 /V ^ a m («) 2; 

^& n (+3V)z _2n ~e "1=1 


(B6) 


n—0 


The coefficients b n (n,N) are polynomials in k and N 
given by 


b n (K,N)= J2 N r '+ r *+- 

ri+ 2 r 2 +.. .-\-nr n =n 

ai(K) ri a 2 (ft) r2 ... a„(>t) r " 
n\r 2 \ ... r n \ 


■ +r n 


(B7) 


where the sum extends over all sets of non negative in¬ 
tegers ri, r 2 , ■ ■. ,r n such that ri + 2 r 2 + ■ ■ ■ + nr n = n. 
For the first five coefficients we get 

bo(K, N) = 1, 

&i(k, N) = | N K (l — 3 K + 2 K 2 ) , 

62 (k, N) = || Nk ( 5 Nn (l — 3 k + 2k 2 ) 2 + 

+6 (-1 + 10 k 2 - 15 n 3 + 6 k 4 )) 
b 3 {n, N) = N k (l — 3 k + 2 k 2 ) x 
x (360 - 18 (-60 + 7N)k + 35N' 2 k 2 - 
- 30 (72 - 42 N + 7 N 2 ) k 3 + 

+5 (216 — 378 iV + 91iV 2 ) k 4 - 
-84 N {-9 + 5 N ) k 5 + 140 N 2 k 6 ) , 

64 («, -/V) = 42525 3V K (l — 3 K + 2 ft 2 ) X 
x (-45360 + 36 (-3780 + 221 N) k- 
-252 (60 — 9 N + 5 N 2 ) k 2 + 

+7 (32400 - 8604 N + 540 N 2 + 25 N 3 ) n 3 - 
-315 (-240 + 36 N - 32 N 2 + 5 N 3 ) n 4 + 

+21 (—10800 + 9684 N — 2700 N 2 + 275 iV 3 ) k 5 - 
-63 (-1200 + 3156 N - 1420 N 2 + 175 N 3 ) k 6 + 

+6 TV (9468 — 10080 N + 1925 N 2 ) k 7 - 
-1260 N 2 (—12 + 5N) k 8 + 1400 N 3 ft 9 ) . 

(B 8 ) 


where 


Now, replacing eq. flB5|) in eq. ( |B2| ) we get 
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1 — /3 r i 1 , K 'P 


-EE 4 " (2k 


- 1 ) 


iV = l n—0 


r (-T) 

x (3 n ^2k - 1 + E) b n (K, N ) ( _ A) -^-D-2n-i = 


OO OO 

= EE^(^) 

JV=1 n=0 


(B9) 


Finally, eqs. (|Bl]) and ( |B9|) lead to the asymptotic ex¬ 
pansion for /'(A)//(A) in eq. (phi). 
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